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A semi-quantitative estimate of the prehadron production time based on recent preliminary HER- 
MES data on hadron transverse momentum broadening in nuclear DIS is presented. The obtained 
production time can well explain the data except for their dependence on the photon virtuality, 
which remains a challenge to current theoretical models. A few mechanisms that may contribute 
to its explanation are suggested, along with possible experimental tests. The average time scale at 
HERMES is found to be of the order of 4-5 fm, comparable to heavy nuclei radii, so that by suitable 
kinematic cuts it will be possible to study both the case in which hadronization starts inside and 
the case in which it starts outside the nucleus. A comparison with recent CLAS preliminary data, 
which qualitatively but not quantitatively agree with HERMES, is performed. 



I. INTRODUCTION 

Nuclear modifications of hadron production in high- 
energy collisions have been observed in both Deep In- 
elastic lepton-nucleus Scattering (nDIS) [1] and in heavy 
ion collisions [2-5]. One typically observes: (i) a suppres- 
sion of hadron multiplicities, called hadron quenching or 
jet quenching; (ii) hadron transverse momentum (pr) 
broadening; (iii) the related modification of the hadron 
PT-spectrum also known as Cronin effect. The nuclear 
modifications can be attributed to the interactions of the 
scattered partons and of the hadrons formed in their frag- 
mentation with the surrounding medium. Experimen- 
tally, partonic in-medium interactions can be isolated 
by studying Drell-Yan (DY) lepton pair production in 
hadron- nucleus collisions. In nDIS and hadron- nucleus 
collisions, the medium is the nuclear target itself, also 
called "cold nuclear matter". In nucleus-nucleus colli- 
sions, the fragmenting parton must also traverse the hot 
and dense medium created in the collisions, be it a hadron 
gas at low energy, or a Quark-Gluon Plasma (QGP) at 
high energy. This medium is also called "hot nuclear 
matter" . 

A precise knowledge of parton propagation and 
hadronization mechanisms obtained from nDIS and DY 
data is essential for testing and calibrating our theoret- 
ical tools, and to determine the properties of the QGP 
produced at the Relativistic Heavy Ion Collider. Con- 
versely, a well known nuclear medium like a target nu- 
cleus allows testing the hadronization mechanism and 
color confinement dynamics in nDIS. Knowledge of par- 
tonic in-medium propagation gained from nDIS, can be 
used in DY scatterings to factor out parton energy loss 
and measure the nuclear modifications of parton distri- 
butions in the initial state [6] . Finally, hadron quenching 
is an important source of systematic uncertainty in neu- 
trino oscillation experiments such as MINOS, which use 
nuclear targets and need to reconstruct the event's kine- 
matic from the hadronic final state [7, 8]. 

Understanding and modeling nuclear modifications of 
hadron production requires knowledge of the space-time 



evolution of the hadronization process [9]. However, 
hadronization is a non perturbative process, and its the- 
oretical understanding is still in its infancy: one has to 
resort to phenomenological models to describe its space- 
time evolution [10-16]. Nonetheless, a few features can 
be expected on general grounds. A parton created in 
a high-energy collision can travel as a free particle only 
for a limited time because of color confinement: it has 
to dress-up in a color- field of loosely bound partons, 
which eventually will evolve into the observed hadron 
wave function. While the naked, asymptotically free, 
parton has a negligible inelastic cross section with the 
surrounding medium constituents, the dressed parton is 
likely to develop an inelastic cross section of the order 
of the hadronic one. Hence, the dressed parton will be 
subject to nuclear absorption in a similar way a fully 
formed hadron is. For this reason, it is usually called 
"prehadron" , and denoted by h* . The prehadron may for 
some time be in a colored state, h% , because color neutral- 
ity is only required for the final state hadron. However, 
it is likely to neutralize its color before hadron forma- 
tion, and the colorless prehadron is denoted h®. We can 
therefore identify 3 relevant time scales, see Fig. 1: (1) 
the "prehadron production time" or "quark lifetime" t p , 
at which the dressed quark develops a sizable inelastic 
cross section, (2) the "color neutralization time" t cn , at 
which gluon bremsstrahlung stops, and (3) the "hadron 
formation time" t^, at which the final hadron is formed. 

For practical applications, hadronization is generally 
pictured as a 2 step process in which the prehadron pro- 
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FIG. 1: Sketch of the time evolution of the hadronization 
process with definition of the relevant time scales. A quark 
q created at time in a hard collision turns into a colored 
prehadron h%, which subsequently neutralizes its color, hi, 
and collapses on the wave function of the observed hadron h. 
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duction time and color neutralization time are identified, 
tp = ten ■ This is a somewhat crude approximation of the 
more complex process sketched above, but is adequate 
to the present status of the theoretical and experimental 
investigation and will be used in this paper. 

In summary, the key quantity we need to investigate 
in order to understand this complex dynamics is the 
hadronization time scale, t p . It is the most general infor- 
mation about the space-time evolution of the hadroniza- 
tion process which can be extracted from experimental 
data. It is the purpose of this paper to provide a semi- 
quantitative and model independent estimate of the pre- 
hadron production time from recent preliminary data on 
the pT-broadcning taken by the HERMES experiment in 
lepton-nucleus scatterings at the HERA accelerator [17 
19]. These data show rich features in terms of their de- 
pendence on the kinematic variables, which are confirmed 
by preliminary data taken by the CLAS experiment at 
Jefferson Lab [20-22] , and are not always simple to inter- 
pret theoretically. I will show to what extent a consistent 
picture of the space-time evolution of hadronization can 
be inferred from these data, and how they are beginning 
to expose the limits of current theoretical modeling of 
this process. 

II. PRODUCTION TIME SCALING 

Hadron quenching in nDIS is studied in terms of the 
multiplicity ratio 
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where Nh is the yield of semi-inclusive hadrons in a 
given kinematic bin, and N e the yield of inclusive scat- 
tered leptons in the same (i/,Q 2 )-bin. The kinematic 
variables z, v, pf, Q 2 are the usual DIS invariants, 
namely, the hadron fractional energy, the virtual photon 
energy, the hadron transverse momentum and the lepton 
4-momcntum transfer squared. 

In Ref. [23] it is conjectured that Rm should not de- 
pend on z and v separately but should depend on a com- 
bination of them: 

Rm = Rm[t{z,v)] , (2) 

where the scaling variable r is defined as 

t = C z x {\- z)v . (3) 

The scaling exponent A can be obtained by a best fit 
analysis of data or theoretical computations. Obviously, 
the proportionality constant C cannot be determined by 
the fit. A possible scaling of Rm with Q 2 is not con- 
sidered because of its model dependence, and because 
of the mild dependence of HERMES Rm data on Q 2 . 
As discussed below, the proposed functional form of r, 



Eq. (3), is flexible enough to encompass both absorp- 
tion models, which assume short production times and 
in-medium hadronization, and energy loss models, which 
assume long lived quarks with (t p ) » Ra 1 where Ra is 
the nuclear radius. The 2 classes of models are distin- 
guished by the value of the scaling exponent: a positive 
A ^ is characteristic of absorption models, while a neg- 
ative A < is characteristic of energy loss models. Thus, 
the exponent A extracted from experimental data can 
identify the leading mechanism for hadron suppression 
in nDIS, and distinguish short from long hadronization 
time scales. 

The scaling of Rm is quite natural in the context of 
hadron absorption models [10-15]. Indeed, prehadron 
absorption depends on the in-medium prehadron path 
length, which depends on the prehadron production time 
{tp), as long as (t p ) < Ra- In the Lund string model [24] 
hadronization is modeled by the breaking of the color 
string stretching from the struck parton to the target 
remnant. The production time is 

z w. 

( tp ) = /(*)(!-*)—£ (4) 

Kstr 

where E q is the struck quark energy, and K s tr the string 
tension. At leading order (LO) in the Strong coupling 
constant a s , the partonic subprocess is 7* + q^q and 
one obtains 

E q = v . (5) 

The factor zE q can be understood as a Lorentz boost 
factor. The (1 — z) factor is due to energy conservation: 
a high- z hadron carries away an energy zE q ; the string 
remainder has a small energy e = (1 — z)E q and cannot 
stretch farther than L = e/K str - Thus the string break- 
ing must occur on a time scale proportional to 1 — z. The 
function f(z) is a small deformation of (t p ), which can 
be computed analytically in the standard Lund model 
[10, 13]. The main features of the estimate (4) are 
dictated by kinematics and 4-momentum conservation, 
hence are of general nature. Indeed they can be obtained 
also by perturbative considerations based on the uncer- 
tainty principle, similarly to what is discussed in Ref. [25] , 
or in the perturbative hadronization model of Ref. [11]. 
The production time (4) is well described by the proposed 
scaling variable r with A > 0. E.g., in the Lund model 
A w 0.7 [23]. In energy loss models [26-28], which assume 
{t p ) » Ra, the scaling is less obvious and holds only ap- 
proximately on a theoretical ground. When performing 
the scaling analysis of the full energy loss models, one 
finds in general A < [23]. 

The central result of the analysis of HERMES data at 
.Blab = 27 GeV performed in Ref. [23] is that pion data 
clearly exhibit 

A w 0.5 ^ . (6) 

As discussed, this is a signal of in-medium prehadron for- 
mation, with production times (t p ) = 0(Ra)- Therefore, 
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FIG. 2: Left: p t -broadening at HERMES for different hadron 
types produced from several nuclear targets as a function of 
the atomic number A (from Refs. [17, 18]). Right: Ratio 
(Pt)a/(Pt)d as a function of xb (from Ref. [19]). The inner 
error bars represent the statistical error and the outer ones the 
quadratic sum of the statistical and systematic uncertainties. 



the scaling variable r can be identified with the average 
production time: 



(7) 



The proportionality constant C in Eq. (3), hence the 
magnitude of the production time, will be estimated in 
Section III. A similar analysis has been attempted in 
Ref. [1], but arbitrarily fixing A = 0.35. The resulting 
Rm(t) shows a rough scaling, with the 2 lowest-z data 
point in each data set deviating from the scaling curve 
and aligning in the vertical direction. This breaking of 
the scaling behavior is more probably due to the par- 
ticular choice of A, rather than to rescattering effects as 
argued by the authors. 



III. pt-BROADENING AND PREHADRON 
FORMATION 

The scaling analysis just described gives only indirect 
evidence for a short production time (t p ), and cannot 
measure its absolute scale. An observable which is more 
directly related to the prehadron production time, and 
allows an estimate of the coefficient C in Eq. (3), is the 
hadron's transverse momentum broadening in DIS on a 
nuclear target compared to a proton or deuteron target 
[11, 29], 
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When a hadron is observed in the final state, neither the 
quark nor the prehadron from which it originates could 
have had inelastic scatterings. Since the prehadron- 
nucleon elastic cross section is very small compared to 
the quark cross section, the hadron's pp-broadening orig- 
inates dominantly during parton propagation. As shown 
in [30, 31], the quark's momentum broadening A(p^) is 
proportional to the quark path-length in the nucleus. If 
the prehadron production time has the form (3) as ar- 
gued in the last section, and as long as the prehadron is 
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FIG. 3: As Figure 2, but for the p t -broadening as a function 
of z (upper panel), v (middle panel), and Q 2 (lower panel). 
Plots taken from Refs. [17, 18]. 



formed inside the nucleus, we obtain 



A(p£) oc (t p ) oc z u ' b (l - z)v 



(9) 



where the exponent 0.5 is determined by the scaling anal- 
ysis discussed in Section II. Then, a decrease of Apj, with 
increasing z or decreasing v would be a clear signal of in- 
medium prehadron formation: indeed, if the quark were 
traveling through the whole nucleus before prehadron for- 
mation, Ap\ would only depend on the nucleus size and 
not on z or v. In addition to the dependence of A(p^) on 
z and v, which primarily stems from energy conservation 
and the Lorentz boost of the hadron, Ref. [11] argues 
that 



A(p 2 T ) 
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The physics behind this proposed behavior is that a 
quark which is struck by a photon of large virtuality ra- 
diates more intensely than for a lower virtuality: as a 
consequence, it will be able to only travel a shorter way 
before hadronization, hence it will experience less px~ 
broadening. A related observable is the Cronin effect, 
which is likewise expected to decrease with increasing z 
or decreasing v : and possibly with increasing Q 2 [11]. 

The HERMES preliminary data on the p-p-broadening, 
reproduced in Figs. 2 and 3 for the reader's convenience, 
show a linear increase with A 1 / 3 , and a clear decrease of 
A(j4) at large z consistent with Eq. (9) and HERMES 
data on the Cronin effect [1]. However, HERMES data 
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TABLE I: Average HERMES kinematics for the p T - 
broadening results [17, 18]. In parenthesis, beside the tar- 
get nucleus symbol is the average in-medium path length 
of the hadronizing system (La) ~ (3/4) .Ra, with Ra = 
(1.12 tmjA 1 ^ 3 . The production time is computed according 
to Eq. (12). The average (xb) is very well approximated by 
(Q 2 )/(m N {u)) [32]. 



also show a clearly increasing A(p 2 -) with Q 2 = 2 — 10 
GeV 2 , at variance with the Q 2 dependence of the pre- 
hadron production time (10) obtained in Ref. [11], and 
with string models in general. They also show an in- 
crease of (pt) \A/(pj,) \D with xb tending to saturate at 
xb > 0.2 [19]. The flat, or even slightly decreasing de- 
pendence of A(p^} on v is also, at first sight, not eas- 
ily interpreted in the context of Eq. (9) and the simple 
hadronization picture used so far. 

For a more quantitative interpretation of the data we 
have first to specify the average kinematics in each bin, 
see Table I and Ref. [18]. Indeed, different observables 
and different experimental bins may be related to differ- 
ent prehadron production times. To this purpose, let us 
consider the A and z distributions: 

• The scaling of A(p^} cx A 1 / 3 up to the Xe nucleus 
indicates that the quark path length is larger than 
the average in-medium path length (Lx c ) of the 
hadronizing system in the Xe nucleon at (z) w 0.4, 
(y) w 14 GeV, so that the prehadron is always 
formed on the surface or outside the nucleus. 

• The above kinematics is close to that of the first 2 
z-bins at z — 0.30 and 0.53. Hence we may assume 
that also in these 2 bins the prehadron is formed 
on the surface or outside the nucleus. However, 
given the decrease of A(p^) with z as z > 0.5, the 
prehadron must soon enter the nucleus at larger z. 
Hence, the production time at z « 0.4 fm cannot 
be much larger than (Lxe)- 

These 2 remarks allow setting the scale for the production 




FIG. 4: z dependence of the prehadron production time 
(t p )/is = tv extracted from HERMES data. 



time in Eq. (3): 
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with z = 0.4 and v = 14 GeV. The average in-medium 
path-length of the hadronizing system can be approxi- 
mated as (L Xe ) « {3/A)R A with R A = (1.12 fm)^ 1 / 3 , 
as if nucleons were uniformly distributed in the nucleus. 
The resulting production time, 



(t p ) « 0.8 z U 5 (l - z) v fm/GeV 



(121 



with v measured in GeV, is plotted in Fig. 4. In principle, 
C should be allowed to depend on Q 2 and Xb, but one can 
neglect it in first approximation at least for the discussion 
of the A-, z-, where (Q 2 ) and (xb) are rather constant. 
With Eq. (12), we can compute the production times for 
each experimental bin, see Table I. 

The ^-distribution is also consistent with this estimate: 
given the production times in Table I, it is not too sur- 
prising that the ^-dependence is basically flat, because 
the prehadron is formed on average on the surface or out- 
side all the studied nuclei. Three effects can contribute to 
tilt the curve and producing the slightly decreasing data: 
(i) a possible dependence of the prehadron cross section 
with v, as it happens for the hadron cross section; (ii) 
medium modifications of the DGLAP parton shower [33- 
36], which can modify the linear dependence of A(p^) on 
the energy loss AE found in [30], and implicit in Eq. (9); 
(iii) the correlation between the py-broadening and the 
average Bjorken variable, (xb), whose origin will be dis- 
cussed in detail in the next Section, and which may in 
fact be the dominant effect. 



IV. pr-BROADENING VS. Q 2 AND x B 

We are now in the position to address the intriguing 
Q 2 distribution. The linear increase of A(p^) with Q 2 
is in sharp contrast with the inverse power dependence 
obtained in the color dipole model, see Eq. (10). On the 
other hand also a flat Q 2 , as would be obtained in most 



■5 



string-based models, seems at variance with these data. 
To my knowledge, no model predicts a substantial rise 
of (t p ) with virtuality: how can we explain these data? 
From the production times shown in Table I, we can see 
that the prehadron is always formed on the surface or 
outside the nucleus, so that the observed Q 2 -dependcncc 
must predominantly have a partonic origin. It cannot 
simply come from multiple parton elastic scatterings, be- 
cause it would be in first order proportional to the in- 
medium parton path-length, which is basically fixed. The 
following 3 mechanisms may contribute to explain these 
data: 

1. Medium- enhanced DGLAP evolution. The DGLAP 
evolution in the kinematics under consideration 
happens entirely inside the nucleus, because the 
prehadron is produced at its surface. The scale Qq 
at which hadronization takes place is not expected 
to strongly depend on the presence of a medium 
[44]. Therefore, in the large-Q 2 bins, a longer and 
medium-enhanced DGLAP evolution [33-36] would 
imply a larger p^-broadening than at low Q 2 . How 
strong this effect is, and whether it can lead to a 
linear increase of the pr-broadening with Q 2 re- 
mains to be seen. This mechanism can be investi- 
gated with a v > v m in cut in the experimental data: 
as f m ; n is increased, the DGLAP evolution would 
increasingly happen outside the nucleus, and the 
slope in Q 2 should become smaller and smaller. 

2. Next-to-leading order processes. The basic hard 
partonic process considered so far is the leading or- 
der in a s quark elastic scattering 7* + q^q. It 
is the dominant process in the valence region at 
large Xb- As xb decreases, however, the gluon dis- 
tribution in the nucleon quickly rises, so that the 
photon-gluon fusion process 7* + + q, which 
is a next-to-leading order (NLO) process, becomes 
non negligible. In this case the photon energy v 
is shared by the quark and the antiquark, so that 
E q < v, therefore reducing the quark lifetime and 
its pT-broadening compared to the LO case where 
E q = v. As can be seen from Table I, the aver- 
age (xb) spanned in the Q 2 distribution lies in the 
transition region between sea partons and valence 
quark dominance. Hence, the competition between 
LO and NLO process might alter the naive picture 
of px-broadening adopted so far, and lead to its in- 
crease with Q 2 . This mechanism may also be at 
work in ^-distributions, although with smaller ef- 
fects because of the more limited range spanned 
by (xb), and contribute to tilt their slope down- 
wards. Note that if this mechanism is indeed at 
work, the estimate of C presented in Eq. (11) only 
gives a lower limit on the hadron production time. 
Experimentally, this mechanism can be tested by 
measuring Q 2 -distributions with suitable cuts on 
xb- The slope in Q 2 should become flat at very 
small or large {xb), see the discussion below. 



3. Colored prehadrons with short formation times. A 
more radical possibility is that at these values of 
(u), the quark does not propagate freely for a long 
time; instead, shortly after the hard interaction, it 
turns into a colored prehadron h% which can loose 
energy by gluon bremsstrahlung, thereby broaden- 
ing its transverse momentum. The prehadron may 
have an inelastic cross section cx 1/Q 2 , growing in 
time to the full hadronic one. If the time evolution 
is slow enough, at low Q 2 the prehadron would be 
subject to more nuclear absorption than at large 
Q 2 , explaining the linear rise of A(p^} with Q 2 
[45] . An extreme version of this mechanism would 
involve 



(*„,) » 0.8 z°- 5 (l - z) v fm/GeV . 

Note that this scenario does not contradict the 
dipolc model prediction that (t p ) cx l/Q 2 , but 
would require the prehadron to propagate as a col- 
ored state for a time of the order of the nuclear ra- 
dius. An experimental investigation of this scenario 
needs very large v, outside the reach of the HER- 
MES experiment but attainable at the Electron- 
Ion Collider (EIC) [37, 38], in order to significantly 
boost (t p ) and allow the quark to propagate as a 
free particle inside the nucleus. 

In order to study the physics behind the experimental 
correlation of A(p^) with xb and Q 2 , one needs to simul- 
taneously address Q 2 -distributions binned in xb and Xb- 
distributions binned in Q 2 . In this way, one can factor 
out the trivial kinematic correlation v = Q 2 j(2m^XB) 
which affects the production time. For example, if 
(tp) oc vjn with k fixed, as in the Lund string model, 
for the LO 7* + q^q scattering we should expect 

Q~ 2 A{p^) XB . hins w const. 
xbA(Pt)q 2 -bins ~ const. 

The deviation from this combined scaling is going to ex- 
pose the underlying dynamics, see Table II: if mechanism 
number 2 is dominant, we would observe Q~ 2 A(p^) XB 
constant in Q 2 but xbA(p\)qi increasing with xb', mech- 
anisms number 1 or 3 would produce an increasing 
Q~ 2 A(pt) Xb but a constant xbA(p^)q2. In the color 
dipolc model [11], t p cx v/Q 2 oc 1/xb hence Q~ 2 A(p^) XB 
would be decreasing and xbA(p^)q2 constant. 

It is also important to extend as much as possible 
the range in over which the p^-broadening is mea- 
sured: at small xb < 0.01 the NLO photon-gluon fusion 
is the dominant partonic process, while at xb > 0.3 — 0.4 
the photon dominantly scatters at LO valence quarks. 
Hence, with mechanism number 2, one may expect A(p^) 
to flat in Q 2 in these 2 regions. 
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TABLE II: Expected results of the scaling analysis of the 
PT-broadening for the models and mechanisms discussed in 
this paper. The numbers in parenthesis refer to the mecha- 
nisms listed in Section IV. The up and down pointing arrows 
indicate increasing or decreasing pT-broadening, the double 
horizontal arrow an approximately constant behavior. 

V. COMPARISON TO CLAS PRELIMINARY 
DATA 

A last question needs to be addressed: how do the 
HERMES data and these scenarios compare to pre- 
liminary CLAS data [20-22] on the PT-broadening and 
Cronin effect? 

At CLAS, the higher beam luminosity allows multi- 
dimensional binning which is only partly accessible at 
HERMES. The kinematics is a bit different due to the 
lower beam energy, -Ebcam = 5 GeV compared to 27 
GeV. The main difference is that (i/) =2 — 5 GeV is 
much smaller than at HERMES, therefore prehadron pro- 
duction should typically happen on shorter time scales, 
mainly inside the nucleus according to the estimate (9). 
Note also that Q 2 = 1 - 4 GeV 2 , and that, typically, 
(xb) clas > {xb) hermes- The main features of CLAS 
data on C, Fe, and Pb targets are the following: 

• The pT-broadening is linear in A 1 / 3 at low A, but 
tends to saturate for large nuclei. The flattening 
is confirmed by the lack of increase in the Cronin 
effect from the Fe target to the Pb target. This 
shows the prehadron being formed on a time scale 
t P £ (4/3)-R Fe = 4.3 fm, smaller than at HERMES. 

• The pT-broadening is rising and then saturating 
with v. The saturation appears at v > 4 GeV, 
where the HERMES px is also approximately flat. 

• The Cronin effect decreases with z as it happens at 
HERMES. 

• The Cronin effect markedly increases as in- 
creases from 0.1 to 0.5. This confirms the Xb de- 
pendence of HERMES pT-broadening. 

All these features corroborate the discussed HERMES 
data. However, the production time I have estimated is 
about a factor 5 smaller than the production time ex- 
tracted in Ref. [29] from CLAS ^-distributions using the 
color dipole formalism [31]. On the other hand, based 
on this paper's analysis, and taking into account the un- 
certainty in the estimate of (La), HERMES data can 
only support up to a doubling of t p , which would oth- 
erwise become incompatible with the large- z decrease of 



the PT-broadening. The HERMES and CLAS data sets 
may be qualitatively reconciled by taking into account 
the observed increase of the PT-broadening with (xb) 
discussed in the previous Section. Indeed, (xb)clas > 
(^b)hermes, which implies a larger pT-broadening at 
CLAS. 

The origin of the discrepancy between the production 
times extracted from HERMES and CLAS data on the 
PT-broadening remains to be clarified. On the theoretical 
side, the two discussed methods to determine the produc- 
tion time should be applied to both data sets, in order to 
check their consistency and further explore the issue. On 
the experimental side, it is very important to establish 
and study the correlation of A(p T ) with xb and Q 2 , as 
emphasized in the previous section. 

VI. CONCLUSIONS 

In this paper, I attempted a semi-quantitative esti- 
mate of the prehadron production time, based on recent 
preliminary HERMES data on hadron pT-broadening 
in nuclear DIS, see Eq. (12). The obtained produc- 
tion time can well explain the dependence of the pr~ 
broadening on A, z and v. Its linear increase with Q 2 
remains a challenge to current theoretical models, and I 
proposed a few mechanisms that may contribute to its 
explanation, also suggesting how to experimentally vali- 
date them. Furthermore, I discussed how the preliminary 
CLAS data qualitatively support the features observed 
at HERMES. 

With the obtained estimate of t p , the prehadrons at 
HERMES are typically formed around the nuclear sur- 
face or slightly outside, while at CLAS they are typically 
formed inside the target. However, a quantitative analy- 
sis indicate that t p \ hermes ~ 0.2£ p |clas- The observed 
x b dependence of the pT-broadening can at least in part 
reconcile the 2 data sets, as well as explain the Q 2 de- 
pendence of the data. It also points at a non-negligible 
role of NLO processes in the hadronization process. 

The measurement of hadron pr-broadening, and the 
related Cronin effect, are beginning to test the limits 
of current theory models on the space-time evolution of 
hadronization, which are based on LO 7*-quark scatter- 
ing, followed by color neutralization of the quark and 
hadron formation. The main theoretical challenges raised 
by these data are, in my opinion, 

• the inclusion of NLO processes in the modeling of 
hadron quenching. 

• the implementation at LO, and subsequently at 
NLO, of alternative models of energy loss, like a 
medium-modified DGLAP evolution; 

• testing of alternative space-time pictures beyond 
the 2 time scale models currently accepted; for ex- 
ample, the role of a colored prehadron should be 
explored. 
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On the experimental side, the proposed physics mech- 
anisms can be verified by measuring the p^-broadcning 
and the Cronin effect with suitable kinematic cuts. In 
particular, 

• the role of parton energy loss can be highlighted by 
large- v/ small- z cuts, such that the parton lifetime 
exceeds the nuclear size; 

• in-medium hadronization can be selected by small- 
^/large-z cuts; 

• the Xb dependence of A(py), the Cronin effect, and 
the multiplicity ratio needs to be better studied 
over a large interval of x B , including small < 
0.01 and large xb > 0.4; additionally, dihadron 
correlations and the hadron multiplicity per event 
as a function of xb can more directly reveal the 
role of NLO processes, in which 2 partons can be 
produced at the hard scattering; 

• the combined analysis of p^-broadening's re- 
distributions binned in Q 2 and of Q 2 -distributions 
binned in is likely to usefully expose the under- 
lying dynamics. 
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